Failure of adaptive self-organized criticality during epileptic seizure attacks 
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Critical dynamics are assumed to be an attractive mode for normal brain functioning as informa- 
tion processing and computational capabilities are found to be optimized there. Recent experimental 
observations of neuronal activity patterns following power-law distributions, a hallmark of systems 
at a critical state, have led to the hypothesis that human brain dynamics could be poised at a phase 
transition between ordered and disordered activity A so far unresolved question concerns the med- 
ical significance of critical brain activity and how it relates to pathological conditions. Using data 
from invasive electroencephalogram recordings from humans we show that during epileptic seizure 
attacks neuronal activity patterns deviate from the normally observed power-law distribution char- 
acterizing critical dynamics. The comparison of these observations to results from a computational 
model exhibiting self-organized criticality (SOC) based on adaptive networks allows further insights 
into the underlying dynamics. Together these results suggest that brain dynamics deviates from 
criticality during seizures caused by the failure of adaptive SOC. 



I. INTRODUCTION 

In the terminology of physics, a system is said to be 
in a critical state if it is poised on a threshold where 
the emergent macroscopic behavior changes qualitatively. 
The hypothesis that the brain is operating in such a criti- 
cal state is attractive because criticality is known to bring 
about optimal information processing and computational 
capabilities [IHJ- Recent experimental observations of 
patterns of neuronal activity exhibiting scale-free dis- 
tributions, a typical hallmark of phase transitions, pro- 
vided further evidence for this hypothesis. Bursts of neu- 
ronal activity were first shown in reduced preparations in 
rat brains to follow power-law probability distributions, 
termed neuronal avalanches [1, 5 . More recently, neu- 
ronal avalanches were also observed in invasive recordings 
from monkeys and cats, strongly suggesting that critical- 
ity is a generic property of cortical network activity in 
vivo [6^-^. 

Additional evidence for the existence of a critical state 
in human brain dynamics comes from a recent contribu- 
tion by Kitzbichler et al. [9 . Using magnetoencephalog- 
raphy (MEG) and functional magnetic resonance imaging 
(fMRI), the authors found power-law probability distri- 
butions of two measures of phase synchronization in brain 
networks. As confirmed by computational models, these 
distributions show power-law scaling specifically when 
those model systems are in a critical state resulting in 
strong evidence that human brain functional systems ex- 
ist in an endogenous state of dynamical criticality at the 
transition between an ordered and a disordered phase. 

Theory predicts local events to percolate through the 
system in the form of avalanches of activity at the criti- 
cal state [10 . Such a critical state requires a homeostatic 
regulation of activity leading to a balance of excitation 
and inhibition in order to prevent states where events are 



either small and local or very large, engaging most of the 
network. A promising mechanism showing robust self- 
organized criticality (SOC) - the ability of systems to 
self-tune their operating parameters to the critical state 
- came from the discovery of network-based mechanisms, 
which were first reported in [11 and explained in detail 
in [121 US] • These works showed that adaptive networks, 
i.e., networks that combine topological evolution of the 
network with dynamics in the network nodes [14 , can ex- 
hibit highly robust SOC based on simple local rules. In 
computational models based on this adaptive interplay 
between network activity and topology, it could be shown 
that local, realistic synaptic mechanisms are sufficient to 
self-organize neuron networks to a critical state, provid- 
ing a plausible explanation of how criticality in the brain 
can be achieved and sustained ^iSHTTj . In summary, ex- 
perimental findings in line with conceptual and compu- 
tational models provide a plausible explanation of how 
criticality in the brain can be achieved and sustained. 

A so far unresolved question concerns the medical rel- 
evance of critical brain activity. Diseases in the central 
nervous system are often associated with altered brain 
dynamics. It has been hypothesized that the dynamical 
properties characterizing a critical state may be seen as 
an important marker of brain well-being in both health 
and disease [18]. Epilepsy is a malfunction of the brain 
associated with abnormal synchronized firing of neurons 
during a seizure [19 . The increased collective neuronal 
firing during attacks has been speculated to be linked to 
a pathological deviation away from a critical state [20] , 
Evidence supporting this idea comes from recent in vitro 
studies of animal brains. There, application of recep- 
tor blockers could drive network dynamics away from its 
normal state where activity patterns of neuron dynamics 
deviated from a power-law ^ |21J . 

Here, we confirm the previously observed critical dy- 
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namics with a complementary experimental methodol- 
ogy, providing additional evidence for the criticality hy- 
pothesis. Furthermore, we show that human brain net- 
works in vivo are not in a critical state during epileptic 
seizure attacks. Deriving the distribution of phase-lock 
intervals (PLI) as an indicator of critical brain dynamics 
from electrocorticogram (ECoG) data, we find the distri- 
bution's scale invariance to be disturbed during seizures. 
Combined with results from a computational model ex- 
hibiting SOC these observations suggest the dynamical 
changes during seizure attacks to rely on a deviation from 
a critical state and hint to the failure of adaptive SOC as 
a cause for seizure generation. Based on the distribution 
of PLI we identify a measure which could prove useful 
for future seizure prediction algorithms. 

II. MATERIALS AND METHODS 

Acquisition And Preprocessing Of Experimental 
Data 

Eight patients undergoing surgical treatment for in- 
tractable epilepsy participated in the study. Patients 
underwent a craniotomy for subdural placement of elec- 
trode grids and strips followed by continuous video and 
ECoG monitoring to localize epileptogenic zones. Solely 
clinical considerations determined the placement of elec- 
trodes and the duration of monitoring. All patients pro- 
vided informed consent. ECoG signals were recorded by 
the clinical EEC system (epas 128, Natus Medical Incor- 
porated) and bandpass filtered between 0.53 Hz and 70 
Hz. Data were continuously sampled at a frequency of 
200 Hz (patients 1-7) and 256 Hz (patient 8, [22 ). The 
study protocols were approved by the Ethics Committee 
of the Technical University Dresden. 
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being a less noisy estimate of Cij averaged over a 
brief period of time At = 2^8 Intervals of phase- 
locking can then be identified as periods when \A(j)ij{t)\ 
is smaller than some arbitrary threshold which we set to 
7r/4 here. 



Estimation Of Deviation Prom A Power-law 
Distribution 

To quantify the deviation from a power-law we defined 
a measure Ap similar to ref. [17]. Ap measures the dif- 
ference between the cumulative density distribution of 
phase-lock intervals and a theoretical power-law distri- 
bution p^^^^ obtained from a fit of the experimental data 
[24] . p^^^^ is calculated from the first time-interval (0- 
150 seconds) of a data set. For each time-interval of 150 
seconds duration, p^^^^ is then subtracted from the cu- 
mulative density distribution of PLI, p^^*^, for each data 
point corresponding to a phase-lock interval x and nor- 
malized by the number of data points N: 

^p = j^^pr'-pT''- (4) 

X 

Positive values of Ap indicate a deviation with in- 
creased intervals of phase-locking, negative values indi- 
cate decreased phase-locking compared to the reference 
power-law distribution. 



Estimation Of Phase Synchronization 

To derive a scale-dependent estimate of the phase dif- 
ference between two time series, we follow the approach 
described in ref. [9 using Hilbert transform derived pairs 
of wavelet coefficients [23 . We define the instantaneous 
complex phase vector for two signals Fi and Fj as: 



CiAt) = 
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(1) 



where Wk denotes the k-th scale of a Hilbert wavelet 
transform and ^ its complex conjugate. A local mean 
phase difference in the frequency interval defined by the 
k-th wavelet scale is then given by 



A^i^j{t) = Arg{Ci^j), 



(2) 



with 



Computational Model 

An infiuential model explaining how dynamical sys- 
tems can self-organize towards a critical state was intro- 
duced in ref. [12 . The mechanism is based on the adap- 
tive interplay between the dynamics of the nodes in the 
network (dynamics on the network) and the rewiring of 
the network's topology (dynamics of the network). More 
precisely, the topology of the network is changed accord- 
ing to the activity of the nodes in the network so that on 
average active nodes loose links and frozen nodes grow 
links. This local rewiring leads to a robust evolution to- 
wards a critical connectivity Kc where the system is at a 
phase transition between order and disorder [12 . 

We first instantiated this model in a network of 200 
randomly interconnected binary elements with states 

= ±1 which are updated in parallel and scanned for 
local rewiring of connections. Under the algorithm de- 
scribed in ref. [12 the network's topology evolves to- 
wards a critical connectivity of Kc ~ 3.1. Our objective 
was to get an estimate of phase-lock intervals between 
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the activity of pairs of nodes for different average con- 
nectivities to provide a reference for comparable analysis 
neurophysiological time series. We therefore monitored 
states (Ji{t) of 20 randomly chosen nodes in random net- 
works instantiated at three different average connectivi- 
ties: K = 2.75, K = 3.1 (critical) and K = 5.0. Each 
network run was limited to 40000 time steps with no 
topological rewiring applied. From these data on the ac- 
tivity in the network cumulative distribution of PLI were 
calculated as described above. We found that the prob- 
ability distribution of phase-lock intervals demonstrated 
power-scaling specifically when the system was at the 
critical connectivity K = 3.1 whereas distributions at 
K = 2.75 /K = 5.0 deviated from a power-law showing 
periods of increased/decreased phase- locking. 



III. RESULTS 

We investigated data sets from ECoG acquired dur- 
ing presurgical monitoring of patients suffering from fo- 
cal epilepsy. Data were continuously sampled at 200 Hz 
(patients 1-7) or 256 Hz (patient 8) with the number 
of channels ranging from 30 to 45 for different patients. 
The time series recorded from the anatomical site where 
the epileptic focus was assumed typically included one or 
more neurographically-identifiable seizure attacks. 

To test brain dynamics for signatures of criticality we 
analyzed ECoG activity in different time windows. The 
data sets were split in intervals of 150 seconds length 
(30000 sample steps at 200 Hz sampling, 38400 in the 
case of 256 Hz) with consecutive intervals overlapping 
by 100 seconds (20000 sample steps at 200 Hz, 25600 at 
256 Hz). Following the approach in [9 , we determined 
the distribution of phase-locking intervals (PLI) as an 
experimentally accessible indicator of critical brain ac- 
tivity. The length of time windows was chosen to be 
long enough to give a reliable estimate of the distribu- 
tion of PLI on the one hand and allow observation of its 
evolution in time on the other hand. For each of these 
sets, we calculated phase- lock intervals and determined 
their cumulative density distributions for scales 2, 3 and 
4 corresponding to frequency intervals 50-25 Hz, 25-12.5 
Hz and 12-6 Hz for patients 1-7 (P1-P7) and 64-32 Hz, 
32-16 Hz, 16-8 Hz for patient 8. 

The distributions for all scales closely follow a power- 
law probability distribution with p{PLI) ~ PLI^ dur- 
ing pre-ictal time intervals. Statistical tests based on 
the Kolmogorov-Smirnov statistic and likelihood ratios 
|24 showed that the hypothesis of a power-law PLI dis- 
tribution could not be rejected for most pre-ictal data 
sets, furthermore a recent comprehensive analysis of var- 
ious fitting functions applied to PLI distributions had 
revealed a power-law to be the most likely fit [9 . The ap- 
parent robustness of the power-law against exact condi- 
tions (different anatomical regions with varying number 
of channels) strengthens the hypothesis of the relevance 
of a critical state in human brain dynamics. 
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FIG. 1. The distribution of phase-locking intervals deviates 
from a power-law during epileptic seizure attacks. Top: The 
electrocorticogram (ECoG) recording shows the onset of a fo- 
cal epileptic seizure attack around 300 seconds time. Bottom: 
Cumulative distributions of phase-locking intervals (PLI) are 
obtained during three time intervals of 150 seconds: pre-ictal 
(left), ictal (middle) and post-ictal (right). Dashed lines in- 
dicate a power- law in the double logarithmic plot. While the 
distribution appears to follow a power-law during the pre- 
ictal period, intervals of increased phase-locking disturb this 
characteristic distribution with the onset of seizure activity. 
Data shown are from patient 1 at scale 3, corresponding to 
the frequency band 25-12.5 Hz. 



While the PLI distribution followed a power-law in 
time intervals preceding the seizure onset, a deviation 
from power-law behavior was observed in intervals con- 
taining the seizure attack. Figure [T] shows distributions 
of PLI derived from a pre-ictal, an ictal and a post- 
ictal time interval. The probability to find longer PLI 
tended to be increased during attacks thereby destroying 
the scale-free property of the original distribution. This 
qualitative change away from a power-law distribution 
could be observed in all 8 patients and across scales (Fig. 

A more quantitative estimate of the deviation from 
the pre-ictal state can be obtained by calculating Ap, a 
measure previously proposed to characterize the diver- 
gence from a critical state [17 . This measure captures 
the deviation from a given empirical distribution from a 
power-law. The power-law fitted to the first (pre-ictal) 
interval was thereby taken as a reference and subtracted 
from the cumulative PLI distributions of subsequent time 
intervals. 

During time intervals preceding the seizure Ap stayed 
at low values indicating no significant deviation from a 
power-law. In time windows containing seizure activity, 
Ap increased to positive values, which is in agreement 
with the qualitative assessment from visual inspection 
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FIG. 2. Comparison of PLI distributions derived from the 
first (pre-ictal) time interval (blue curve) and an interval dur- 
ing the seizure attack (red curve). Distributions from seizure 
intervals tend to exhibit longer periods of phase-locking re- 
sulting in a deviation from a power-law of the distribution's 
tail. Plots are shown for scale 3 corresponding to the fre- 
quency band 25-12.5 Hz for patients 1-7 (P1-P7) and 32-16 
Hz for patient 8 (P8), respectively. 



showing a divergence from the initial distribution. After 
seizure attacks, a slow decrease of Ap could be observed 
suggestive of a relaxation process back toward a power- 
law distribution (Fig. [3|. 

Although Ap values hinted that the power-law distri- 
bution was conserved during pre-ictal intervals, we ob- 
served variations in the amount of phase-locking as time 
approached the seizure onset. To investigate this with a 
better resolution in time we repeated the calculation of 
PLI for shorter time windows of 5000 sampling steps (cor- 
responding to 25 seconds for patients 1-7 and 19.53 sec- 
onds for patient 8) every 50 sampling steps (correspond- 



ing to 0.25 or 0.19 seconds). While the power-law be- 
havior of the PLI distribution appeared to be conserved 
until seizure onset, we found that the number npLi and 
sum of all phase- lock intervals ^ PLI varied during pre- 
ictal periods. Figure [4] shows the product of the two, 
npLi ■ ^PLI^ as a time series for each patient. Fur- 
thermore, this measure tended to decrease toward the 
time of seizure onset. For the scales 2-5 under investi- 
gation, the decrease was occurring most prominently in 
the frequency band of scale 4 (12-6 Hz P1-P7 and 16-8 
Hz P8). While the amount of phase- locking varied and 
showed the tendency to decrease toward seizure onset 
times, the power-law tail of the PLI density distribution 
was conserved indicating to a robust underlying mecha- 
nism producing the power-law distribution. 

For obtaining further insights into the underlying dy- 
namics of the power-law probability distribution of PLI 
and its absence during epileptic seizure attacks, we com- 
pared experimental results to a simple computational 
model exhibiting self-organized criticality. The net- 
work model introduced by Bornholdt and Rohlf [12 self- 
organizes toward a critical state and through its simplic- 
ity allows for an understanding of the underlying mech- 
anism by which this self-tuning is accomplished. Specif- 
ically, the adaptive interplay of network dynamics and 
topology, a mechanism also at work in more elaborate 
models of SOC in neural networks p!5HT7] , robustly or- 
ganizes systems parameters, in this case the average con- 
nectivity toward values Kc where the network's state 
is at a phase transition between ordered and disordered 
dynamics. 

At the self-organized connectivity, statistical tests [24] 
revealed that the hypothesis of a power-law for the dis- 
tribution of PLI can be accepted whereas it has to be re- 
jected for distributions at connectivities below or above 
Kc (Fig. [5|. There, the distribution of PLI deviated 
from a power-law consistent with a state away from crit- 
ical dynamics. The distribution at connectivities corre- 
sponding to the ordered phase of network dynamics is 
shifted towards larger PLI similar to the one observed 
during epileptic seizure attacks (bottom right in Fig. [5|. 

The close agreement between patient and model data 
suggests that the deviation from a power-law observed 
during epileptic seizure attacks relies on a shift of dy- 
namics toward an ordered phase. It further hints that it 
is the mechanism of adaptive SOC, the ability to tune 
system parameters to values where network dynamics is 
at a phase transition and PLI are distributed according 
to a power-law, that fails during epileptic seizure attacks 
in neuron networks in the brain. 



IV. DISCUSSION 

The relevance of critical brain dynamics is currently a 
heavily debated topic. Indirect evidence for such a state 
comes from power-law distributed observables in neuro- 
physiological data. Although various mechanisms can re- 
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FIG. 3. Development of the deviation from a power- law. ECoG recordings from 8 patients showing a focal seizure attack are 
shown along with Ap values for consecutive time windows of 150 seconds duration overlapping by 100 seconds. The power-law 
fit of data in the first time window was taken as the reference to calculate Ap. Although different in extent, an increase of Ap 
quantifying the deviation from the initial pre-ictal distribution can be observed during seizures for all patients and different 
scales. 



suit in an event size distributions exhibiting power-laws 
[25] , such distributions arise when a system is in a critical 
state [TO] . Recently, the power-law distribution of phase- 
lock intervals between pairs of neurophysiological time se- 
ries was shown to be a specific hallmark of dynamic crit- 
icality in human brain dynamics [9 . Using this indicator 
on ECoG data, a complementary experimental method- 
ology to [9], we confirm the previously observed critical 
dynamics, providing additional evidence for the critical- 
ity hypothesis. Secondly, we show that the critical state 
is disturbed during epileptic seizure attacks. More pre- 
cisely, the distribution of the PLI synchronization mea- 
sure deviates from a power-law, characterizing the criti- 



cal state of normal neuronal dynamics, during epileptic 
seizures, providing the first direct evidence of disturbed 
critical dynamics related to a pathology in vivo. 

Our findings support the notion of a physiological de- 
fault state of balanced brain dynamics between regimes 
of exuberant and frozen activity. Physiological neuronal 
activity is characterized by intermittent periods of syn- 
chronization between different anatomical regions. In 
terms of dynamical system's theory, such a state cor- 
responds a critical state at a phase transition between 
order and disorder. A deviation from this balanced state 
toward dynamics with pathologically increased times of 
synchronous activity as observed in epileptic patients 
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FIG. 4. Evolution of the amount of phase- locking quantified 
by ripLi • ^PLI, with upli being the number and ^PLI 
the sum of all phase-lock intervals. The measure varies dur- 
ing the time course and tends to decrease reaching a mini- 
mum (marked by black arrows) shortly before the beginning 
of the seizure. Seizure onset times are marked by vertical 
dashed lines. Phase- locking intervals were determined from 
overlapping time windows of 5000 sampling steps length for 
the frequency band defined by scale 4. 
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FIG. 5. Distribution of PLI in a model exhibiting self- 
organized criticality. Through an adaptive interplay of net- 
work dynamics and topology, the Bornholdt model evolves 
toward a critical state with an average connectivity Kc- The 
left graph shows the evolution to Kc ^3.1 for a network with 
200 nodes for different initial connectivities. Right: cumula- 
tive distributions of PLI for scale 1 in random networks at 
average connectivities K = 5.0 (green), K = 3.1 (blue) and 
K — 2.75 (red). The PLI distribution follows a power-law at 
the self-organized connectivity Kc = 3.1. 



leads to a deviance from the physiological critical state 
resulting in impaired functionality. 

Optimal functional capabilities of neuron networks 
have been related to a critical state before [Tl[26l|27]. In 
this work we observed that power-law scaling of phase- 
lock intervals is conserved until the seizure onset at which 
the power-law character of the PLI distribution is lost rel- 
atively abruptly. This is remarkable since other measures 
characterizing the PLI distribution such as the sum and 
number of all PLI are varying and tending to decrease 
prior to the seizure. With power-law scaling as an in- 
dicator of criticality this suggests that a state of critical 
dynamics is maintained until the seizure onset and lost 
thereafter which manifests in impaired functionality and 
the start of clinical symptoms. 

In comparison to the relatively late loss of the charac- 
teristic power-law distribution, other precursors are much 



earlier indicative of a seizure onset. The identification of 
precursors of a seizure attack is a question of high clin- 
ical relevance [28 . Earlier studies reported a decrease 
of synchronization-based measures suggesting them as 
promising indicators of a pre-ictal state [29-32J. We also 
observed a tendency of decreased phase-locking prior to 
the seizure attack. In a recent study using a measure 
based on phase-locking the authors found pre-ictal de- 
creases of synchronization levels to be occurring mostly 
within the 4-15 Hz band [31 similar to our findings which 
were most evident for scale 4. Given the amount of in- 
formation on phase-locking in a specific frequency band 
contained in our measure ripLj - ^PLI one can specu- 
late that it should perform better in detecting pre-ictal 
changes of neural dynamics than other less specific syn- 
chronization measures. However, it should be noted that 
in order to assess the usefulness and performance a sta- 
tistical validation taking into account the sensitivity and 
specificity of this measure is needed which is beyond the 
scope of this article. 

A mechanism by which complex networks can self- 
organize toward a critical state is based on the adaptive 
interplay between the dynamics on the network, i.e. neu- 
ronal activity, and the dynamics of the network, i.e. the 
shaping of synaptic connections. Through this interac- 
tion system parameters can be locally tuned to a state 
of global criticality [121 US] • While the simple model de- 
scribed in this work captures these essential ingredients 
allowing for an understanding of the underlying concept, 
more elaborate mechanism can be expected to be at work 
in real- world neuron networks [TSHIZl ES] • It is conceiv- 
able that physiological neuron networks in the brain tune 
their parameters to more than one parameter to reach a 
state of criticality. Besides the average connectivity of 
the network K., the balance between excitation and inhi- 
bition, for example, has been shown to be an important 
parameter to sustain a homeostatically balanced critical 
state and prevent regimes of overly synchronized activ- 
ity. It is conceivable that it is the robust mechanism of 
adaptive SOC maintaining neuron networks in the brain 
close to a critical state characterized by dynamics ex- 
hibiting power-law probability distributions even while 
network dynamics undergoes changes as captured by the 
variations of overall phase-locking, for example. 

The deviation from a power-law distribution of PLI re- 
ported here corresponds to a shift away from a balanced 
critical state and to our knowledge constitutes the first 
direct proof of impaired critical dynamics related to a 
pathology in vivo. This observation is supported by ex- 
perimental results from in vitro studies. The application 
of receptor blockers in slice preparation of animal brains 
resulting in an excess of excitation in the network de- 
stroyed the power-law distributed avalanches of neuronal 
activity and led to increased avalanche sizes correspond- 
ing to a super-critical state O [6] . Analogously, human 
tissue removed from epilepsy patients exhibited abnor- 
mally regulated avalanches with periods of hyperactivity 
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In summary, experimental results from in vitro experi- 
ments [5j E] and in vivo observations presented here com- 
bined with insights from computational models based on 
adaptive SOC fTH [TSHTTI l33] suggest the failure of the 
adaptive interplay between neuron activity and network 
topology to lead to the deviation from a critical state 
with pathological, in the case of epilepsy overly synchro- 
nized, activity patterns. The sudden loss of the power- 
law distribution during the onset of an epileptic seizure 
attack hints to some threshold after which adaptive SOC 
is not able to maintain a critical state and fails. A devi- 
ation from the default critical state towards a dynamical 



regime with decreased phase-locking is also conceivable. 
For instance in neurodegenerative diseases with impaired 
neuronal connectivity, the deviation from a power-law of 
PLI could potentially be used to identify and characterize 
pathological conditions. 
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